Centrality Dependence of Direct Photons in Au+Au Collisions at ^snn = 200 GeV 
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We calculate the spectra of high energy photons emitted in relativistic Au+Au collisions for 
various centralities and compare to data recently collected at the Relativistic Heavy Ion Collider by 
the PHENIX collaboration. Our results for photons from primary hard scatterings and photons from 
interactions of jets with the medium are consistent with the measurements of neutral pion and direct 
photon production in p + p collisions and give a good description of direct photon spectra measured 
in Au+Au collisions. The contribution of photons from jet-to-photon conversion in the medium can 
be as large as the photon yield from hard scatterings in the momentum range px « 2 ... 6 GeV/c. 
We show that this novel mechanism is not ruled out by any existing data. 



High energy nuclear collisions are studied to under- 
stand the properties of strongly interacting matter at the 
highest energy densities and to search for the transition 
of hadronic matter to a quark-gluon plasma . Electro- 
magnetic probes are excellent tools to gather information 
about the dense and hot medium created in such colli- 
sions. The mean free path of photons is much larger than 
the typical size of the emerging fireball, allowing them to 
escape nearly unperturbed once created in the system 
0|- Since photons are emitted throughout the history 
of the fireball, one expects to obtain information that is 
complementary to that gained from the measurement of 
hadrons, which are subject to strong final-state interac- 
tions and are reflecting the latest stage of the collision. 
Thermal radiation of photons can occur both from a 
uark-gluon plasma 0, 0| and from a hot hadronic gas 
Because the photon spectrum reflects the thermal 
distribution of the emitting matter, it can serve as a ther- 
mometer. On the other hand, high energy photons are 
most likely to be created either in hard processes between 
partons in the colliding nuclei , or in final-state interac- 
tions of high energy jets 0. Photons from jets naturally 
carry information about the medium they traverse and 
hence can serve as a medium probe. 

The latter process was introduced in our earlier work 
where we predicted the existence of a medium induced 
component in the photon spectrum in the intermediate 
Pt range. In this novel mechanism, a hard quark pass- 
ing through the medium converts into a photon carrying 
most of its transverse momentum. Recently Zakharov 
has developed the related concept of modified photon 
bremsstrahlung of jets in the medium j^. An up-to-date 
calculation of photon emission including such jet-medium 
interactions can be found in the work of Turbide et al. 
|lO|. The general arguments given for photons in this 
work also apply to the emission of lepton pairs initiated 
by quark jets in the plasma 

Relativistic nuclear collisions at different centralities 
result in the formation of systems of different sizes, initial 
temperatures and life times. Assuming that the yield of 



processes involving large momentum transfer (hard pro- 
cesses) scale as the number of collisions A^coii, while those 
with small momentum transfer (soft processes) scale as 
the number of participants iVpart, one can predict the 
relative importance of the two contributions as a func- 
tion of centrality (see e.g. 0|). The recently published 
data on single photon production in p + p collisions 
at y^snn = 200 GeV [13 and the centrality dependence 
of the photon yield in Au+Au collisions measured by 
PHENIX at the Relativistic Heavy Ion Collider (RHIC) 
[1^ is our first chance to obtain an understanding of the 
different photon sources. 

Here, we shall address the question whether the data 
allow us to draw conclusions about the relative contri- 
butions from different photon production mechanisms to 
the direct photon spectrum in Au+Au collisions. The 
data is still limited by (statistical and systematic) uncer- 
tainties, but the results represent a remarkable achieve- 
ment because of the experimental challenge to subtract 
the background of photons from secondary hadronic de- 
cays, especially neutral pions. Fortunately, in Au+Au 
collisions at RHIC this background is reduced at inter- 
mediate and high transverse momenta p-r > 2 GeV/c by 
the large final state suppression (commonly called "jet 
quenching") for pions and other hadrons [lajllal- Thus, 
while the low-pT part of the photon spectrum remains 
clouded by the large hadronic background, this is an in- 
vitation to look for electromagnetic signals of a quark 
gluon plasma at intermediate and high pr 

Let us briefly recall the photon sources that are im- 
portant at > 2 GeV/c. We define prompt photons as 
those coming from hard Compton and annihilation pro- 
cesses (e.g. q+g (7+7 and q+q — > 5+7) and those emit- 
ted from jets after the hard interaction as bremsstrahlung 
(e.g. q ^ q + j)- A basehne for these processes can be 
obtained by studying p+p and p+p interactions. Besides 
an obvious scaling from p + p collisions to Au+Au col- 
lisions with the number A^coib and an isospin correction 
due to the presence of p + n and n + n collisions, Comp- 
ton and annihilation spectra are expected to be altered 
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only marginally by initial state nuclear effects. Neutral 
pion spectra measured in d+Au collisions indicate that 
such effects from shadowing and fcT^-broadening are small 
[l^ . The hard process itself is pointlike on typical scales 
of the medium and therefore remains unaffected. Direct 
photon spectra measured in d+Au could help to quantify 
the role of initial state effects. 

Bremsstrahlung emission from jets, on the other hand, 
takes place on more extended distances and thus can 
be affected by final state interactions. Energy loss 
effects would then be noticeable as a suppression of 
bremsstrahlung [Tflj ]. while secondary scattering in the 
medium could initiate additional bremsstrahlung or mod- 
ify bremsstrahlung by coherence effects 0. In our work 
we assume that the total photon yield from jets can 
be described by the superposition of the usual vacuum 
bremsstrahlung with an additional contribution from jets 
annihilating and Compton scattering in the medium as 
described in Q. 

Photons from secondary hard scatterings in the pre- 
equilibrium phase have been studied elsewhere and could 
contribute to photons at intermediate and high pT |2Cll| . 
In order to minimize possible effects from neglecting pre- 
equilibrium emission, we assume a very early thermaliza- 
tion time, tq ~ 0.15 fm/c, in central collisions, not much 
larger than the formation time of minijets. An early ther- 
malization of the matter created in Au+Au collisions at 
RHIC is supported by several observations, especially the 
large elliptic flow. We also assume that the matter is 
chemically equilibrated at the initial time. Secondary 
collisions are then either between thermalized partons or 
between hard scattered partons and the medium, which 
are both taken into account here. Similar strategies for 
dealing with the uncertainties of p re-equilibrium emission 
have been invoked elsewhere |2 ij . 

We now recall the specifics of the jet-medium inter- 
action process. We denote the momenta of the jet, the 
thermal parton, and the photon by Pjot, Pth and p^, 
respectively. The leading order QCD Compton and an- 
nihilation cross sections are peaked in the forward and 
backward directions. In the laboratory frame we have 
|pjet| ^ |pth| ^ T, where T is the temperature of the 
plasma. For high energy photons, i.e. |p^| ^ T, this im- 
plies that |p-y| « |Pjct|- This justifies calling the process a 
conversion of a jet (or more precisely: a quark initiating 
a jet) into a photon with similar momentum 8J. Only 
quark jets contribute to this process, because the gluon- 
photon scattering cross section is not peaked at forward 
angles. 

The rate of photon production by annihilation and 
Compton scattering of jets in the medium can be ap- 
proximated as [3 
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where C = -1.916, = 5^T^/6 and = 5^/(47r) is 
the strong and a the electromagnetic coupling, q denotes 
all light quark and antiquark species with charge , and 
fq is the distribution of minijets of flavor q. It is worth 
emphasizing that the conversion property of the process 
is reflected in ^ by the fact that the photon spectrum is 
directly proportional to the jet spectrum fq. Calculations 
without some of the the simplifying approximations Q 
that lead to the expression in Eq. iQ produce only small 
modifications of the photon yield |1Q|] . 

In order to calculate the yield from thermal emission 
and jet conversion, we need to model the space-time evo- 
lution of the fireball. We want to compare with data 
taken at midrapidity, hence we assume the boost invari- 
ant Bjorken scenario for the longitudinal expansion. In 
addition we neglect the marginal transverse expansion of 
the plasma on the jet-photon conversion part, i.e. during 
the short time jets are propagating through the QGP. 
But the transverse expansion is fully taken into account 
for the thermal contribution. 

We integrate over the conversion vertex x in (|T)) in 
terms of the space-time rapidity ry, the radial coordinate 
r (with r = |r|) at which the jet is emitted, and the 
proper time r, so that d^x = rdrdrid^r. We parameterize 
momenta p by their rapidity y and transverse momentum 
Pt, {pt = |pt|). Starting from the yield dNjd/ (d^prdy) 
of quarks initiating a jet, the distribution of hard partons 
in a boost-invariant scenario is given by 0, |22| 



fqiP,x) 



(27r)3 



JCt 



gqirR^TpT d'^prdy 



p{r)6{y-ri) 



(2) 



where Rt is the radius of the fireball, gq = 6 the spin- 
color degeneracy factor of parton q, and p is a transverse 
profile function normalized to one. For simplicity, we will 
take the profile to be azimuthally symmetric with area 
density p{r) = 2(1 — r'^ / R^)Q{Rt — r) and transverse 
radius Rt = l.2{N^^,t/2Y'\ 

The number of participants, iVpart, is a function of cen- 
trality, but we neglect the azimuthal asymmetry in non- 
central collisions. This approximation does not allow us 
to predict the azimuthal distribution of photons, but it 
introduces only a small error for the yield integrated over 
the azimuthal angle. The jet spectrum for an arbitrary 
impact parameter h is obtained by scaling 
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This expression shows that the photon yield due to jet- 
photon conversion scales as a product of the nuclear 
thickness and a term dependent on the initial conditions, 
like temperature and size, and on the history of the evolu- 
tion of the fireball. A measurement of the photon spectra 
at different centralities could thus help reveal this infor- 
mation, as the primary hard photons will scale with the 
nuclear thickness Taa alone. 
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FIG. 1: Invariant cross section Eda/dPp for tt" and 7 pro- 
duction in p + p collisions at — 200 GeV as a function 
of transverse momentum pr at v — 0. Data points are from 
the PHENIX experiment [lllgg. We also show a full NLO 
calculation |27|. 



We use a leading order perturbative calculation of jets 
and of direct Compton and annihilation photons , us- 
ing CTEQ5L parton distributions |2^. For gold nuclei 
EKS98 corrections of the nuclear parton distributions 
[2^ are taken into account. Higher order perturbative ef- 
fects are accounted for by a phenomenological i^-factor, 
which is known to depend weakly on the transverse mo- 
mentum Pt at RHIC energies [23 and is here assumed 
to be constant. 

Vacuum photon bremsstrahlung is obtained by convo- 
luting the jet distributions with the photon fragmenta- 
tion functions given in [3 (A = 200 MeV), and we use 
KKP fragmentation functions ^261 for the yield of neutral 
pions. We calculate the 7r° and photon cross sections at 
midrapidity for p+p and compare to the values measured 
by PHENIX. This allows us to fix the K factors for jet 
production, -ftTjet, and direct Compton and annihilation 
photons, K^, independently. We obtain a good descrip- 
tion of the p + p data with Kjct — 1-8 and = 1.5, see 
Fig. These values agree with those reported in ref . 
on the basis of a next-to-leading order (NLO) perturba- 
tive QCD calculation and our results also match the NLO 
calculation by Gordon and Vogelsang . 

Once the if-factors are fixed we apply them to cal- 
culate the yield of direct photons in Au-f Au collisions. 
The prompt photon yield in Au-f Au is close to that in 
p + p scaled with Taa , due to small shadowing correc- 
tions. Note, however, that there is an effect from the 
admixture of p + n and n + n collisions with different 
photon yields. Therefore a simple scaling of p+p results, 
as often invoked by the experimental collaborations for 
comparison, introduces a small but systematic bias of the 
photon yield from hard processes towards larger values. 

Next we use the jet yields in Au-|-Au to calculate the 
photons from jet-photon conversions using and |(2Jl. 



Centrality 


0-10% 


10-20% 


20-30% 


30-40% 


40-50% 


Npmt 


326.2 


234.5 


166.0 


114.0 


75.0 


Taa [1/mb] 


22.75 


14.35 


8.00 


5.23 


2.86 


TO [fm/c] 


0.151 


0.166 


0.182 


0.201 


0.221 


To [MeV] 


0.434 


0.395 


0.361 


0.328 


0.297 



TABLE I: Number of participants A'part, nuclear thickness 
factor Taa, initial time tq and initial temperature To for dif- 
ferent centrality bins. 



For the fireball we assume a thermally and chemically 
equilibrated plasma at some early time tq with initial 
temperature Tq. The product T^t is conserved during 
the isentropic longitudinal expansion, and its value can 
be fixed by the total hadron multiplicity dNh/dy. The 
initial conditions are determined by imposing the ther- 
malization condition tq w 1/(3To). Tabled lists the ini- 
tial times and temperatures inferred from the charged 
hadron multiplicity dN/drj measured by PHENIX (1^ 
for different centrality bins. Nuclear thickness factors 
and number of participants are also given. 

The transverse profile of the initial temperature is 
fixed by assuming that the energy density scales with the 
square of the nuclear thickness factors, so that Ti{r) = 

To [2(1 - rVi?T)] This has been found to work well 
for the study of the centrality dependence of hadronic 
spectra [2^ as well as lepton pair production 30], and 
should be sufficient for our purpose. Starting from the 
initial time tq, we integrate to either the time t^, 
when the medium has cooled to the critical temperature 
Tf — 160 MeV, or to the time when the jet reaches the 
boundary of the fireball, whichever comes first. 

For the calculation of thermal photons we assume that 
the thermally and chemically equilibrated quark-gluon 
plasma, which undergoes a boost-invariant longitudinal 
and an azimuthally symmetric transverse expansion 31], 
converts to a hadronic gas below Tc- The production of 
photons from the deconfined phase is calculated using 
the complete leadiirg order emission rate 4]. The latest 
results of Turbide et al. for photon emission from 
hadronic matter are used. 

Figs. El El El and El present our results for photon yields 
at midrapidity in Au-|-Au collisions as a function of trans- 
verse momentum in four different centrality bins: 0-10%, 
10-20%, 20-30% and 40-50%. Both prompt photons and 
photons from jet-plasma interactions are shown, and the 
sum is compared to PHENIX data ,14]. For completeness 
we also show our result for thermal photons. The agree- 
ment with data is generally quite good for all centrality 
bins. This is also true for the other centrality bins, up to 
the 60-70% bin, which are not shown here. 

As already pointed out in 8], the jet-photon conver- 
sion spectrum is falling off with px faster than the spec- 
trum of prompt photons, a consequence of the additional 
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FIG. 2: Photon yield dN/ {cPpTcLy) as a function of pt for y = 
in central (0-10%) Au+Au collisions at ^ = 200 GeV. We 
show primary hard photons (dashed), jet-photon conversion 
(solid and labelled) and the sum of both (upper-most solid 
curve). Thermal photons are also shown (solid and labelled). 
Data are from the PHENIX collaboration 0]. 



FIG. 5: The same as Fig.Hbut for the 40-50% centrality bin. 
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FIG. 3: The same as Fig.Elbut for the 10-20% centrality bin. 
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FIG. 4: The same as Fig.Hbut for the 20-30% centrality bin. 



hard momentum transfer "with a propagator ~ l/pj'. In- 
terestingly the jet-photon conversion yield is as large as 
both the primary hard photon yield and the thermal yield 
around pr = 2 GeV/c in central Au-|-Au collisions. Un- 
fortunately, the quality of data in this region is still poor. 
Bet"ween 4 and 6 GeV/c, "where the data provide a better 
constraint, jet-medium photons still add roughly a 50% 
contribution on top of primary hard photons. 

It has been argued that photons from primary hard 
processes alone are sufficient to describe the Au-|-Au data 
[l^ Elf and that no additional contribution is needed. 
The first statement is correct, to a certain extent, as can 
be seen in Fig.|2| Ho"wever, its value is limited by the size 
of the experimental error bars. We have demonstrated 
here that a calculation taking into account jet-photon 
conversion in a consistent manner also leads to a result 
"which is compatible "with present data. Caution should 
be excercised -when isospin effects are omitted in the dis- 
cussion. 

Energy loss effects for jets before they convert into pho- 
tons have been investigated by Turbide et al. 10]. The 
net effect on the yield of photons from jet-medium inter- 
actions is found to be small, about 20%. This is because 
most photons are emitted at a time "when the plasma 
is still hot. Jets have then traveled only a short dis- 
tance through the plasma and have not lost a significant 
amount of energy. 

The centrality dependence of the data is "well described 
by our calculations. While direct photons scale "with the 
number of collisions A^coib photons from jet-medium in- 
teractions exhibit a slightly stronger scaling. The relative 
importance of this process thus decreases in more periph- 
eral collisions. This can also be seen in Fig. "where we 
plot the nuclear modification factor 



Raa = 



dNAA/dy\p^ye GcV/c 
iVcoU dNpp/dy\p^yQ gcV/c 



(4) 



calculated from the photon spectra integrated for px > G 
GeV/c vs centrality. We sho"w our calculation and a ver- 
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FIG. 6: The nuclear modification factor Raa for direct pho- 
tons and pT > 6 GeV/c as defined in I^J as a function of 
the number of participants A^part (solid line) . We also show a 
calculation effectively taking into account energy loss of jets 
before convert into photons (dashed line). Data are taken 
from PHENIX collaboration Q. 



sion taking into account an effective 20% energy loss of 
jets before conversion into photons and compare with 
PHENIX data [l4j. Due to the large error bars the data 
is compatible with 1 for central and midperipheral colli- 
sions, in accordance with a binary collision scaling from 
p + p collisions. Our calculation is compatible with the 
data as well, and predicts a slight rise of Raa with cen- 
trality. Our results show that the importance of the con- 
tribution from jet-photon conversion to the total photon 
spectrum grows below px — G GeV/c. A high statistics 
measurement of Raa in this region would be very useful. 

In summary, we have calculated the photon spectrum 
in Au-|-Au collisions at RHIC energies resulting from pri- 
mary hard photons, jets interacting with the medium and 
thermal radiation. The calculations are consistent with 
neutral pion and photon spectra measured in p + p. We 
obtain a good description of the px dependence and cen- 
trality dependence of photon production in Au-|- Au. The 
contribution from jet-photon conversions can be as large 
as 100% (50%) that of primary hard photons at 2 (5) 
GeV/c in accordance with data measured by PHENIX. 

At the Large Hadron Collider (LHC), the importance 
of the jet-photon conversion process will be significantly 
enhanced. It has been shown that it will be the dominant 
contribution below pT ~ 12 GeV/c [Rlioj. 

This work was supported by DOE grants DE-FG02- 
87ER40328 and DE-FG02-05ER41367. 
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